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We have measured the dependence of the relative integral cross section of the reaction 
Li+HF— >-LiF+H on the collision energy using crossed molecular beams. By varying the intersection 
angle of the beams from 37° to 90° we covered the energy range 25meV< En < 131 meV. We 
observe a monotonous rise of the cross section with decreasing energy over the entire energy range 
indicating that a possible translational energy threshold to the reaction is significantly smaller than 
25 meV. The steep rise is quantitatively recovered by a Langevin-type excitation function based on 
a vanishing threshold and a mean interaction potential energy oc R~^'^ where R is the distance be- 
tween the reactants. To date all threshold energies deduced from ab-imtio potentials and zero-point 
vibrational energies are at variance with our results, however, our findings support recent quantum 
scattering calculations that predict significant product formation at collision energies far below these 
theoretical thresholds. 



I. INTRODUCTION 

The rapid development of techniques for cooling, trap- 
ping and manipulating cold and ultracold molecules 
opens the opportunity of studying chemical reactivity in 
the low temperature regime [IHl]. Cold reactive colli- 
sions are not only important limitations to stable trap- 
ping of molecules but also present fascinating aspects of 
their own interest. In the cold collision regime, reactions 
are governed by quantum dynamics involving tunneling 
and scattering resonances [21 [3 • For our study of cold 
collision phenomena we have selected the title reaction 
due to both the availability of theoretical results and the 
applicability of a variety of experimental methods. 
Attractive from the theoretical point of view is the small 
number of electrons together with three chemically very 
different but light atomic constituents that make the 
slightly exoergic (wO.lGmeV) reaction Li+HF ^-LiF+H 
an ideal prototype system for developing methods to 
calculate reliable ab-initio potential energy surfaces 
(PESs) [6 16^ and to calculate quantum scattering phe- 
nomena jS] [12 dni fT7H52] . Particularly challenging for a 
quantitative prediction of observables is the low energy 
range near threshold where quantum effects dominate 
and a precise knowledge of the PES around the transition 
state is essential. Crucial for the assessment of the pre- 
ciseness and reliability of the applied theoretical methods 
is the comparison of experimental and computational re- 
sults. Very qualified for this purpose are certainly angu- 
lar and velocity distributions of products |15[ [33] but also 
the dependence of the relative integral reaction cross sec- 
tion (IRCS) on the translational collision energy, Etr, the 
excitation function [221 [23 [MIIMl [SS] . At low collision 
energies the excitation function bears direct information 
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on the existence and size of a translational threshold that 
allows a sensitive examination of the shape of the PES 
near the transitions state and the height of a possible 
potential energy barrier. 

An inspection of computational results shows that all 
PESs |6HT6| known to date feature a significant bar- 
rier at a bent transition state with substantial height. 
However including zero point energies of both the reac- 
tants and the transition state the semi-classical trans- 
lational threshold energies deduced from these barrier 
heights turn out to be markedly smaller. For example, for 
the most frequently used recent PESs, the semi-classical 
thresholds and barriers amount to «27 and 182 meV |12j . 
« 68meV and 233 meV [TS], «56 and 221 meV [13 [16], 
respectively. Neglecting quantum scattering phenomena 
the semi-classical thresholds manifest the lowest colli- 
sion energies required for product formation from ground 
state reactants. An observation of products at ener- 
gies below these calculated thresholds could be a clue to 
marked quantum phenomena or to insufficient accuracy 
of the PESs or to a combination of both. Recently, it has 
been demonstrated in time-dependent wave packet scat- 
tering calculations [5: 32 based on the PES of ref. [15 , 16] 
that resonances play an important role for Li+HF. They 
lead to product formation at a significant rate at collision 
energies far below the semi-classical threshold. 
One interesting experimental aspect of the title reaction 
is the availability of methods that allow to investigate 
scattering phenomena in a wide range of collision en- 
ergies from hot (several 100 meV) to cold (20meV) and 
ultra-cold (< 1 meV). In an earlier experimental study we 
have measured the excitation function within the range 
82meV < Etr < 376 meV using molecular beams 

intersecting at a 90° angle. The energy was varied by 
applying the seeded beam technique. At energies above 
« 120 meV the function turns out to be roughly constant 
but below it rises steeply with decreasing Etr and assumes 
at the lowest energies the largest values within the entire 
range. Consequently, a threshold to reaction, if existing 
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at all, has to be located below 82 meV. 
In the present paper we report the results of another 
crossed beam study on the excitation function designed 
to lower the range of translational energy to 25 meV. Vari- 
ation of the collision energy is achieved by varying the in- 
tersection angle of the beams (see below). We find that 
in the extended energy range the IRCS continues to rise 
monotonously with decreasing £'tr thus shifting the up- 
per boundary of a threshold to below 25 meV. 
It should be noted that a novel experimental setup us- 
ing new techniques for preparing samples of cold atoms 
and molecules is about to be completed. It will allow to 
continue our search for the threshold of this process at 
energies below 1 meV. 



II. EXPERIMENTAL METHODS 




FIG. 1: Schematic view of the molecular beams arrangement. 
The setup is basically the same as used in our previous exper- 
iments |35| except for the HF beam that can now be rotated 
around the intersection volume to vary the intersection angle 
F in the range 37° < F < 90°. An additional HF-beam mod- 
ulator placed between nozzle and skimmer is not shown for 
clarity. 

To achieve translational energies below 82 meV (the low- 
est energy accessible in ref. |35] ) we abandon the usual 
90° geometry and cross the beams at a variable intersec- 
tion angle [36"38] . Fig.[l]shows a sketch of the molecular 
beams arrangement. It deviates from the one used in ear- 
lier experimental studies |551[551I55] only by the setup of 
the HF-beam. The nozzle is now mounted within a differ- 
entially pumped vacuum chamber (not shown) that can 
be rotated around the scattering volume. In this way the 
intersection angle F of the two beams and thus also the 
mean relative velocity 

V = [vli + <F - ^vuVHF COS F]°-5 (1) 



TABLE L Beam velocity parameters, velocities in m/s 



Run liLi OLi uhf ohf 

A/B 1530 221 1205 185 
C 1750 290 1215 195 



and the translational collision energy 

E^,.iT) = ^^lv' (2) 

can be varied continuously, denotes the reduced mass 
of the colliding reactants and vu, v^f are their nom- 
inal velocity. Depending on F, significantly smaller or 
larger collision energies can be achieved compared to the 
one obtained for the usual perpendicular intersection. In 
this study the two vacuum chambers housing the beam 
sources are designed such that F can be varied in the 
range 37° < F < 90°. 

The experimental method and data acquisition follows 
closely the ones discussed in ref. [35]. Briefly, the neat 
HF beam is created by a nozzle heated to around 500 K 
to avoid clustering. Stagnation pressure (250mbar) and 
nozzle temperature are always kept constant. The ele- 
vated temperature leads to a mean rotational energy of 
8.9 meV [3S] and a population of the first vibrationally 
excited state of < 0.25%. The beam is chopped with a 
frequency of a few Hz using a tuning-fork like modula- 
tor or a chopper wheel at a duty cycle of 50% and the 
intensity is monitored by a quadrupole mass spectrome- 
ter detector equipped with an electron bombardment ion 
source. The Li beam is diluted with Ne and its intensity 
is measured by surface ionization on a hot rhenium (Re) 
ribbon. During an experimental run only the intersection 
angle F is varied while the operational conditions of the 
beams are kept constant. 

The velocity of both beams is determined by conventional 
time-of-flight (TOF) arrangements. The (density) veloc- 
ity distributions, n{v), are extracted from the measured 
TOF profiles by fitting the parameters u, a of the ex- 
pression 

n{v) = const exp{ — [(u — w)/a]^} (3) 

to the data. The parameters u and a of the various data 
sets are compiled in Table 1. The velocity spreads of the 
Li and HF beams a result in an experimental relative 
uncertainty of the collision energy Etr of about 25%. 

The products created in the intersection volume of the 
beams are detected via surface ionization on a hot Re 
ribbon mounted in a separately pumped ultra high vac- 
uum chamber (residual gas pressure below 10~*mbar). 
A channeltron converts the ions desorbing from the Re 
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surface to electron pulses which are counted by a two 
channel scaler synchronized with the HF beam modula- 
tor. One scaler counts the signal and background, the 
other only the background pulses. The detected scat- 
tering intensity (signal) /tot is then derived as difference 
of the two scaler contents. Angular distributions of the 
signal /tot (6) £^re obtained by rotating the main detec- 
tor automatically around the intersection volume in the 
plane of the beams within a wide range of laboratory 
(LAB) scattering angles Q. Crucial for the data analy- 
sis (see below) is the knowledge of time-of-flight (TOF) 
distributions of scattered particles at various scattering 
angles. They are measured employing a fast spinning 
chopper wheel with 8 equally spaced slots (2 mm wide) 
mounted between the skimmer of the HF beam and the 
scattering volume (see Fig. [ij . The length of the flight 
path is 254 mm. 
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FIG. 2: Typical TOF profiles measured at F = 37° and 50° 
for tlie indicated scattering angles. The profiles are used to 
separate elastic (smooth curves) from the detected scattering 
intensity (data points). See text for more details. 

The surface ionization is not specific with respect to alkali 
atoms or alkali compounds and detects both species with 
roughly the same efficiency. Therefore, the detected scat- 
tering intensity /tot is proportional to the sum of fluxes of 
both the elastically scattered Li atoms /lj and the prod- 
ucts LiF, /lif- To separate both components we employ 



the TOF distributions; examples are displayed in Fig. [2] 
for r = 37° and 50°. All TOF proflles feature a peak at 
short flight times and a broad shoulder or even a second 
peak at longer times. At F = 90° and higher Li velocities 
the peaks are always well separated |35j . 



(a) r=37 




FIG. 3: Newton diagrams illustrating the kinematics for three 
beam intersection angles. The radius of the outer circle 
around the tip of the centroid velocity vector C corresponds 
to the velocity of the elastically scattered Li atoms in the CM 
system, the radius of the inner filled circle to the maximal 
product velocity. 

The collision kinematics are illustrated in Newton- 
diagrams (Fig. |3| relating the LAB and center-of-mass 
(CM) frame velocities of reactants and products for beam 
intersection angles F — 37°, 50°, and 90°. The light ar- 
rows represent the laboratory frame velocity vectors vli 
and vhf of the reactants Li and HF which include the in- 
tersection angle F. The radius of the outer circle around 
the tip of the centroid velocity vector C corresponds to 
the velocity of the elastically scattered Li atoms in the 
CM system, the radius of the inner shaded circle repre- 
sents the maximal product velocity. The light straight 
lines that are tangent to the inner circle indicate the 
range of scattering angles AO where reaction products 
LiF are expected to occur in addition to elastically scat- 
tered Li atoms. Elastically scattered Li atoms along C 
have higher speeds than the LiF products and will there- 
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fore reach the detector at shorter flight times. Li atoms 
that are elastically scattered in the opposite direction in 
the CM frame have low speeds in the LAB frame and will 
spatially disperse on their way to the detector such that 
no significant contribution to the signal is expected. 
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FIG. 4: Angular distributions of the detected (upper dia- 
grams) and of the product (lower diagrams) scattering inten- 
sity measured at F = 37° and 50°. The open circles in the 
upper diagrams mark the elastic intensity deduced from the 
TOF distributions of Figure 2. The solid line is the best fit 
interpolation used to deduce the reactive scattering intensity 
(lower diagrams) by subtracting the elastic from the detected 
scattering intensity. 

Thus, in Fig. [2] the peak at short flight times (fast parti- 
cles) is attributed to elastically scattered Li atoms while 
the one at later times (slow particles) is associated with 
reaction products LiF. The solid line represents a numer- 
ically simulated TOF-distribution of elastically scattered 
Li atoms whose reliability was confirmed by a comparison 
with TOF profiles measured for the pure elastic system 
Li+Ne that features practically the same kinematics as 
Li+HF (open circles in Fig. The elastic TOF pro- 
file is then fitted to the fast slope of the data. The area 
below this scaled elastic peak over the area below the en- 
tire TOF profile eventually provides the fraction of the 
detected scattering intensity associated with elastically 
scattered atoms. For a given F this fraction has been 
determined at scattering angles for which TOF distribu- 
tions are available and used to isolate the angular distri- 
bution of the elastic component /Li(0) from the detected 
scattering intensity /tot(0)- The result of this procedure 
is shown in Fig. [4] (upper panel) . The open circles rep- 
resent the angular distribution of the elastic component 
/Li (6) and the solid line is a best fit curve through these 
points. The difference between the detected scattering 
intensity and the elastic component eventually yields the 
angular distribution of the product flux /LiF(0) (lower 



panel in Fig.|4|. For more details see |35| . 
In principle the signal contribution identified with re- 
active scattering could also be generated by rotation- 
ally inelastic collisions. At F 50° inelastic collisions 
Li+HF (j = 0) -^Li+HF(j = 5) would yield Li atoms 
with similar CM speeds as the LiF products. At F = 37°, 
Li atoms could be generated at speeds similar to those of 
LiF in inelastic collisions Li+HF (j' = 0) — >-Li+HF(j' = 
3). Such selective rotation-changing collisions appear 
quite implausible, though, given the common scaling laws 
that predict a fast decay of the rotational inelastic colli- 
sion cross section with increasing level spacing |40| . 



III. RESULTS 

The IRCS, (Tr(Etr), is proportional to the total flux of 
products TVlI'p^' generated in the scattering volume V and 
defined by the expression 



^L?F"V("Li nHF V V). 



(4) 



Here, riLi, ?^HF and v denote the number densities of 
the indicated beams at the intersection volume and the 
mean relative velocity, respectively. Deviating from our 
earlier study |35], we leave the operational conditions of 
the beam sources constant and vary only the intersection 
angle F. Thus both densities njjF and riLi are constant, v 
can be easily deduced from the most probable beam ve- 
locities (eq. [T|) , and an inspection of the intersection ge- 
ometry shows that V = Vgo" / sin F holds approximately. 
The crucial quantity iV^^p^' is not directly accessible in 
the present in-plane scattering experiment but can be 
deduced from the measured total in-plane product fiux 

xin— plane • r i 

^LiF using the lormal expression 



total 
LiF 



7\rtotal 
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^LiF 



(5) 



Inserting eq. |5] into eq. |4] and suppressing all constant 
quantities we obtain the expression 



rin— plane 
'LiF 



sinV/v /l"f 



in— plane 



(6) 



relating the measured quantities and the relative integral 
reaction cross section. The ratio of fluxes in eqs. |5]and[6] 
corrects for the fraction of products that miss the detec- 
tor; it can be readily deduced from the relative differen- 
tial reaction cross section (DRCS) in the center-of-mass 
frame. In a previous study [33J using perpendicularly in- 
tersecting beams we have measured the relative DRCS 
at 6 energies within the range 82meV< Eu- < 376 meV 
and found that the ratio is constant with respect to the 
energy within an error margin of ±4% |35] . Unfortu- 
nately, a comparable extensive investigation of DRCSs is 
not yet available for the present low energy range. How- 
ever, a preliminary analysis of the product angular distri- 
butions measured at F = 37° (24 meV) and 50° (45 meV) 
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indicates that from 119 meV (F = 90°) to 24meV 
(r = 37°) a transition from the forward/backward to a 
preferred sideways type DRCS occurs. Taking this into 
account we find for all data sets that the product of fac- 
tors left to the in-plane flux in eq. [6] is constant within 
the band width ±5%. In view of a forthcoming more so- 
phisticated determination of the DRCS we suppress these 
small corrections and derive the relative IRCS from the 
data using the simplified relation 

ar{E,,) (X (7) 

We have measured three sets of angular distributions A, 
B, and C at various intersection angles for the beam pa- 
rameters given in Table 1. The constancy of the oper- 
ational conditions was checked carefully by measuring a 
reference angular distribution repeatedly during one run. 
The intersection angles range from 37° to 90° (A,B) and 
48° < r < 90° (C) corresponding to the energy range 
25meV < Etr < 108 meV (A,B) and 50meV < Etr- < 
131 meV (C). For each angular distribution of a given 
set we measured between four and six TOF-profiles to 
separate elastic and reactive scattering. 
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FIG. 6: The excitation function deduced from data set A 
(circles), B (triangles) and C (diamonds). The data sets are 
mutually normalized to obtain the best agreement within the 
overlap. The solid line through the data points is given by 
the simple power law (7r(i5tr) oc (l/iStr)"'^. 
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FIG. 5: Product angular distributions of data set A mea- 
sured at the intersection angles 90° (open square), 82° (open 
circle), 75° (open triangle), 68° (open diamond,) 62° (solid 
square), 55° (solid circle), 45° (solid triangle) and 37° (solid 
diamonds) . Note the shift of the curves to smaller angles with 
decreasing F caused by the kinematics (Fig. [3| and the strong 
increase of the peak intensity. The area below the curves is 
the total product intensity used to determine the IRCS. 

As an example we show in Fig. [5] the product angular dis- 
tributions of set A. With decreasing F the curves shift to 
smaller LAB angles according to the changing kinemat- 
ics (Fig. |3| and their peak intensities rise dramatically. 

The total in-plane intensity /l'if^''*"'' '-'^ fr(£'tr) (eq.M) is 
given by the sum over all product intensities multiplied 
by the Lab-angle increment (area below the curves), 

^in^pla„c^^^^.^(Q.)^Q.^ (8) 



The results for set A are compiled in Table 2 together 
with those for B and C and displayed in Fig. [6] as a func- 
tion of the collision energy. The data points are nor- 
malized such that they agree optimally within the over- 
lapping energy range. Their statistical error is not in- 
cluded in Fig. |6] for clarity. It can be estimated from 
the scatter of the points and amounts to about ±5% es- 
sentially due to uncertainties occurring in the process of 
separating elastic and reactive scattering. The solid line 
through the data points represents the simple power law 
ariEt,) oc (l/i?tr)°-^ 



IV. DISCUSSION 

A strong motivation for performing the present experi- 
ments was the search for an answer to the question raised 
in the earlier study |35j : How does the excitation func- 
tion continue below a collision energy of 82meV? The 
previous results indicated two possibilities: either the 
excitation function assumes a maximum followed by a 
decline towards a threshold as one would expect for a 
reaction featuring a non-vanishing translational thresh- 
old energy or it continues to increase as expected for a 
reaction without threshold. The results of the present 
investigation strongly support the assumption that no 
threshold hinders the reaction of Li with HF. The excita- 
tion function continues to ascend monotonously at least 
down to 25meV and is likely to continue this way [41J. 
Due to the elevated nozzle temperature of 500 K, required 
to suppress dimerisation, low rotational states of HF 
are populated with a mean rotational energy of 8.9 meV 
corresponding to a mean rotational quantum number of 
j = 1.4 [35]. Thus the indicated threshold energies refer 
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TABLE II: Total in-plane scattering intensities in arbitrary 
units 
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FIG. 7: Excitation function including the results of our previ- 
ous study |35| (solid squares). The new and earlier results are 
mutually normalized to obtain the best agreement within the 
overlap. The curves at energies < 200 meV are Langevin cross 
sections based on potential energy functions oc (dotted), 
oc (dashed) and oc R~^'^ (solid). The solid smooth line 
for energies > 200 meV refers to a constant hard sphere col- 
lision model matched to the Langevin functions at 200 meV. 
The light solid line represents the result of recent wave packet 
calculations |32) . 



on average to the reaction Li+HF {v = 0j'=1.4). This 
internal energy may increase the threshold of the ground 
state reaction Li+HF(u — 0,j = 0) on average by the 
mean rotational energy. 

In Fig. [7] we have combined the results of the present 
and earlier experiments and obtain an excitation function 
that covers now a wide energy range of more than one 
order of magnitude from 25meV to 376 meV. The curve 
can be subdivided into three sections: above 200 meV 
it is more or less constant, below 130 meV the function 
rises monotonously with decreasing iStr, and in between 
it gradually changes the shape correspondingly. This be- 
havior suggests the following classical reaction model. 

(a) At energies > 200 meV a rigid sphere collision mech- 
anism |34l I42| prevails where the reaction occurs with 
a constant, energy independent probability V whenever 
the impact parameter b is smaller than the maximal one 
bmax- The resulting constant IRCS is then = Virbf-^^^^. 
The mechanism neglects long range forces and thus all 
trajectories of the approaching reagents are straight lines. 

(b) Below 130 meV long range forces start getting im- 
portant. The force curbs the trajectories of the colliding 
reagents toward the center of the force and thus colli- 
sions with impact parameters larger than 6niax may re- 
act. With decreasing energy the influence of the forces 
grows and thus the maximal impact parameter and the 
IRCS increase. Provided there is a vanishing threshold 
the IRCS continues to rise and eventually diverges classi- 
cally for i?tr 0. Assuming that (i) the interaction po- 
tential energy is cx 1 /i?" where R is the distance between 



the reagents, and (ii) the reaction occurs with a constant 
probability whenever the reagents overcome the effective 
potential's maximum, then the IRCS is described by the 
"Langevin" power law (T,{Etr) oc {l/Etr)^/' [51111]. The 
best fit curve in Fig. |6] corresponds accordingly to an in- 
teraction potential oc In Fig.|7]the s — 2.5 result 
(solid line) is compared with curves referring to s = 6 
(dotted line) and s = 4 (dashed line) for energies below 
200 meV. The corresponding potentials describe the long 
range forces of an atom-molecule (s = 6) and of an ion- 
molecule system (s = 4). Both curves feature the typical 
steep rise but are clearly at variance with our data. The 
best fit power s — 2.5 indicates that in our low energy 
range the chemical forces near the hard sphere radius 
are responsible for the rise of the IRCS rather than the 
asymptotic long range forces with s — 6. 
(c) Between 130 and 200 meV the transition between the 
hard sphere and the potential energy dominated mecha- 
nism occurs. 

The prediction of this classical model is illustrated in 
Fig. [7] as solid line where we have omitted the transi- 
tion region and matched directly the hard sphere and 
Langevin functions at 200 meV. The curve recovers the 
data nearly quantitatively over this wide energy range. 
The above classical model suggests a vanishing thresh- 
old but the experimental results provide only an upper 
boundary and a small threshold may well be in agree- 
ment with our findings. According to classical mechanics 
such a translational energy threshold is tightly related 
to the height of the potential energy barrier Vt, separat- 
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ing reagents and products. Provided the approaching 
molecules are in their vibrational and rotational ground 
states {v — 0, j = 0) the formation of products requires 
that the translational energy of the approaching particles 
exceeds Vh- A non vanishing reaction threshold is then 
the consequence of an existing barrier. 
Quantum mechanics relaxes the tight classical relation 
between threshold and barrier height. The colliding par- 
ticles possess zero-point vibrational energy that varies 
during the approach adiabatically from the zero-point 
energy of the free molecule E^^^^ to the one of the tri- 
atomic aggregate at the barrier E^. Taking this into 
account, product formation for ground state molecules 
{v = 0,j = 0) is then allowed in a classical sense, if the 
sum of the collision and zero-point energy of the reagent 
molecule exceeds the sum of the potential and zero-point 
energy at the barrier. The resulting semi-classical thresh- 
old energy 



(9) 



includes in addition to the barrier height Vh also the rel- 
evant zero-point energies E^ and E^^^. However, in con- 
trast to classical mechanics, eq. [9] constitutes no rigorous 
lower limit for the translational energy leading to prod- 
uct formation. Due to quantum phenomena such as tun- 
neling and resonances products may be formed also at 
collision energies far below this semi-classical threshold 
(see below). 

To date all PESs available for Li+HF predict a significant 
barrier height Vh ^-'16 . Including the zero point energy 
of HF (256 meV) and of the transition state (95±10meV 
depending on the geometric structure of the PES) the 
semi-classical threshold (eq. [9| for the most recent PESs 
amounts to 27meV [12J, 56meV [HI [16] and 68meV [T3] . 
The present experiments provide as upper boundary for 
the threshold energy 25 meV and, considering the signif- 
icant product flux observed for the Ar seeded Li exper- 
iment IS], 17meV. Furthermore, the excellent fit of the 
data predicted by a model based on a vanishing thresh- 
old suggests a threshold located at energies ^ 17meV. 
Thus our results suggest that the experimental threshold 
energy lies markedly below the predicted semi-classical 
values. This discrepancy could be due to insufficient ac- 
curacy of the computational methods used to calculate 
the PESs or to quantum effects that promote reactions at 
energies far below the semi-classical threshold. The latter 
is supported by recent wave packet calculations based on 
the PES of |15l [TB] that predict a threshold for the total 
reaction cross section at about lOmeV which is quali- 
tatively reconcilable with our measurements [32J. The 



theoretical excitation function [32] is displayed in Fig. 7] 
The curve exhibits only a slight modulation around the 
value (T,. sa O.5A as the translational energy is reduced 
down from Etr ~ 250 meV to a threshold value of about 
lOmeV (see Fig. [t]). The shape of the predicted excita- 
tion function departs significantly from the experimen- 
tal one but the threshold energy is in accord with the 
present result. As the main source for product forma- 
tion near threshold collisions with small J (total angular 
momentum) have been identified. The reaction proba- 
bilities for J = exhibit a rich spectrum of oscillations 
and resonances as a function of energy [5] [3^ but these 
structures disappear if realistic IRCSs are calculated by 
summing over all J. 



V. CONCLUSION 

The benchmark reaction Li+HF— >LiF+H has been stud- 
ied at translational energies down to 25 meV using a new 
crossed-beam apparatus with variable scattering angle 
between Li and HF beams. This arrangement allows 
to tune the translational energy while keeping the beam 
source conditions constant. The integral reactive scat- 
tering rate (excitation function) is deduced from angle- 
resolved scattering as well as from time-of-fiight traces 
in comparison with purely elastic Li+Ne scattering. The 
resulting excitation function, which extends earlier mea- 
surements to lower energies, steeply rises as the colli- 
sion energy falls below about 150 meV. This behavior 
is consistent with a barrier-less Langevin-type reactive 
process with i?~^-^-scaling of the atom-molecule interac- 
tion potential. Alternatively, an energy threshold below 
^ 25 meV may be present as predicted by recent wave 
packet simulations |32| . 

Future efforts to further reduce the collision energy for 
conclusively disclosing or ruling out the existence of a 
reaction threshold will demand more sophisticated ex- 
perimental approaches. To this end, a magneto-optical 
trap for preparing an ultracold Li scattering target is cur- 
rently being set up and will be combined with a source 
for slow and cold molecules based on a rotating nozzle 
and electrostatic guiding [15] . 
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